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ABSTRACT 



Context. BA-type supergiants are amongst the most optically-bright stars. They are observable in extragalactic environments, hence potential 

accurate distance indicators. 

Aims. An extensive record of emission activity in the Ha line of the BA supergiants jSOrionis (Rigel, B8Ia) and a Cygni (Deneb, A2Ia) is 

indicative of presence of localized time-dependent mass ejections. However, little is known about the spatial distribution of these apparent 

structures. Here, we employ optical interferometry to study the Ha line-formation region in these stellar environments. 

Methods. High spatial- (~0.00r') and spectral- (R=30 000) resolution observations of Ha were obtained with the visible recombiner VEGA 

installed on the CHARA interferometer, using the S1S2 array-baseline (34m). Six independent observations were done on Deneb over the 

years 2008 and 2009, and two on Rigel in 2009. We analyze this dataset with the ID non-LTE radiative-transfer code CMFGEN, and assess 

the impact of the wind on the visible and near-IR interferometric signatures, using both Balmer-line and continuum photons. 

Results. We observe a visibility decrease in Ha for both Rigel and Deneb, suggesting that the line-formation region is extended (~1.5- 

1.75/?*). We observe a significant visibility decrease for Deneb in the Sin637lA line. We witness time variations in the differential phase 

for Deneb, implying an inhomogeneous and unsteady circumstellar environment, while no such variability is seen in differential visibilities. 

Radiative-transfer modeling of Deneb, with allowance for stellar-wind mass loss, accounts fairly well for the observed decrease in the Ha 

visibility. Based on the observed differential visibilities, we estimate that the mass-loss rate of Deneb has changed by less than 5%. 



Key words. Techniques: high angular resolution - Techniques; interferometric - Stars: emission-line - Stars: mass-loss 
(HD 197345, HD 34085) - Stars: circumstellar matter 



Stars: individual 



1. Introduction 

Supergiants of spectral types B and A (BA-type supergiants) 
are evolved massive stars of typical initial mass of 25-40 M© 
and high luminosity (>1O^L0). Their luminosity and tem- 
perature place them among the visually brightest massive 
stars. Therefor e, they are particu larly interesting for extragalac- 
tic astronomy (IPuls etal.l 120081) . Moreover, they represent at- 



tractive distance indicators by means of the identified wind- 
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momentum-luminosity relationship (WMLR) (IKudritzki et al 



20081: lEvans & Howarthl200llKudritzki et al.ll999h . As a con 
sequence, nearby BA supergiants have been analyzed with 
sophisticated radiative-transfer tools. Two subjects of intense 
scrutiny, both observationally and theoretically, have been 
Deneb (g Cygni, HD 197345, A2 Ia) and Rigel (^Orionis 



HD 34085. B 8Ia) (ISchiller & Przvbilla2008clAufdenberg et al 



2006bll2002l for Deneb). 

Intensive spectroscopic monitoring of the activity of wind- 
forming lines, such as Ha, suggests that the stellar-wind vari- 
ability of luminous hot stars is associated with localized and 
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co-rot ating surface structures. The pioneering work of iLucv 
(119761) based on a period analysis of radia l- veloc ity curves 
of Deneb obtained in 1931/32 by ' Paddockl d 19351) revealed 
the simultaneous excitation of multiple non-radial pulsations 
(NRPs). This eff ort was follo wed by the extended optical 
Heros campaigns (iKaufenl 19981) . monitoring for ~100 consec- 
utive nights late B-type/A-type supergian ts (Kaufer et al. 1997 , 



1996allbh and early B-type hypergiants (IRivinius et alJll997h . 



Attempts to associate this activity with large-scale surface 
magnetic fields have so far been unsuccessful ( Schnerr et al. 
2008; Ve rdugo et alj2003l:lHenrichs et alj2003l;lBvchkov et al.1 
2003i) . However, the appearance of thin surface convection 



zones at the supergiant stag e may favor the forrn ation and emer- 
gence of magnetic fields dCantiello et al.ll2009l) , and may also 
have a significant impact on potential NRPs. 

BA-type supergiants such as Rigel (B8 la) and Deneb 
(A2Ia) exhibit observational evidence of the random and 
pseudo-cyclic activity of the stellar wind. A pulsation-driving 
mechanism has often been proposed, although the associ- 
ated micro-variabil ity is obse r ved fo r a wide range of lumi- 
nosity and log Tgff. 'Gautschr ('1992') argued for the presence 
of strange modes. Gautschy (2009) tentatively proposed that 
Deneb's micro-variability is caused by a thin convection zone 
efficiently trapping the non-radial oscillations, and significant 
progress has bee n made in deter mining the mechanism leading 
to mass ejection ( ICranmen2009h . The collective effect of multi- 
ple NRPs was also proposed as a promising means of explain- 
ing the large value of the macro-turbulence parameter; NRPs 
also repre sent an attractive mechanism for the f ormation of Be- 
star disks dAerts et alJl2"009HNeiner et al.ll2009l) . 

Based on intensive monitoring of spectroscopic lines, many 
observational campaigns have attempted to discriminate be- 
tween these models. Monitoring of Deneb and Rigel is still on- 
going (JDa nezis et al. 2009; Morrison et al. 2008; Rzaev 2008; 



Markova et a l. 2008; Markova & Puis 2008i: lMorelet al.l i2004; 



Rivinius et aljil997i ; .Kaufer et aljil996h, for a non-exhaustive 



list). 



The uniform-disk (UD) angular diamet ers of Deneb and 



Rigel ar e 2.4 mas and 2.7 mas, respectively dAufdenberg et al. 
[2008; M ozurkewich et al.l l2003i). making them good targets 



for accurate long-baseline optical interferometry. An exten sive 
study of Deneb was performed bv lAufdenberg et al.ld2002l) . us- 
ing different radiative-transfer models (with allowance for the 
presence of a wind), and constraints from optical interferome- 
try data. Deneb's wind should enhance the limb darkening rel- 
ative to hydrostatic models that neglect it. However, this can 
only be constrained observationally by using long baselines in 
the near-lR (longer than 250m) or by observing in the visible, 
thus relaxing the constraints on the baselines by a factor of 
about 3-4. The Ha line is particularly interesting in that con- 
text. Its large opacity makes it very sensitive to changes at the 
stellar surface and above, in particular through modulations in 
stellar-wind properties. Its extended line-formation region is 
therefor e more easily resolv e d than the deeper-forming con- 



tinuum. iDessart & Chesneaul d2002h investigated the possibil- 
ity of monitoring the mass-loss activity in wind-forming lines 
by means of optical interferometers equipped with spectral de- 
vices of sufficient dispersion. 



BA supergiants rotate slowly (v sin / of about 25-40km s"' ), 
at least relative to their terminal wind-velocity of ~200- 
400km s '. Spectrally resolving the Doppler-broadened 
Balmer lines thus requires a resolving power as high as R =10 
000. At optical wavelengths, the most conspicuous indicator 
of mass-loss in early-type supergiants is Ha, which exhibits 
the distinctive P-Cygni profile morphology and has often been 
a crucial estimator of the mass-loss rate in spite o f its doc- 
umented variability. While lAufdenberg et al.l (120021) was suc- 
cessful in reproducing many observables obtained for Deneb, 
the modeling of the Ha line proved difficult. The most sig- 
nificant and persistent discrepancy between the synthetic and 
observed profiles has been the depth of the absorption com- 
ponent, which is significantly weaker in the observed spec- 
trum (residual intensity ~0.6). Furthermore, while the veloc- 
ity of the absorption component minimum is quite variable, it 
is rarely, if ever, shifted blueward by mo re than ~50km s~\ 
which is ~20% of the terminal velocity. ISchiller & Przvbilla 



(2008) used a different strategy, using the Ha profile as a ref- 
erence for a detailed modeling of the wind characteristics, and 
reached a more satisfactory solution. Nevertheless, they also 
reported some remaining discrepancies in matching the Ha ab- 
sorption, and proposed additional effects that may need to be 
considered in the modeling, such as wind structure or the influ- 
ence of a weak magnetic field. 

The VEGA recombiner of the CHARA array is a recently 
commissioned facility that provides spectrally dispersed inter- 
ferometric observables, with a spectral resolution reaching R 
=30 000, and a spatial resolution of less than one mas. The in- 
strument recombines currently the light from two telescopes, 
but 3-4 telescope recombination modes are foreseen in a near 
future. The Ha line of bright, slow rotators such as Deneb 
or Rigel can be isolated from the continuum, and the spatial 
properties of the line-forming region can thus be studied with 
unprecedented resolution. Using the smallest baseline of the 
CHARA array (baseline of 34m), we conducted a pioneering 
temporal monitoring of Deneb uncovering a high level of ac- 
tivity in the Ha line-forming region. Rigel was also observed a 
few times. 

The paper is structured as follows. In Sect. |2]we present 
the optical interferometry data and their reduction, in addition 
to complementary spectra obtained by amateur astronomers. 
We then review in Sect. 13. li the previous interferometric mea- 
surements aiming to more tightly constrain the diameter of 
Rigel and Deneb. We then apply this extensive record of mea- 
surements in the visible and near-lR continuum to interpret 
semi-quantitatively the spectrally-dispersed measurements in 
Sect. 13. 21 and to elaborate a radiative-transfer model of the stars 
in Sect. H 

2. Observations and data processing 

2.1. VEGA observations 

Deneb and Rigel were observed wit h the Visible spEctroG raph 
And Polarimeter (VEGA instrument, M ourard et al.l2009l) inte- 
grated within the CHARA array at Mount Wils on Observatory 
(California, USA. Iten Brummelaar et al.ll2005h . Deneb obser- 
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Table 1. VEGA/CHARA observing logs. 



Date 


Time 


SCI/CAL 


Projected Baseline 








Length [m] 


PAn 


2008.07.28 


06:04 


Deneb 


33.4 


6.2 


2008.07.28 


09:01 


Deneb 


32.9 


-18.8 


2008.07.28 


10:51 


Deneb 


26.2 


-46.7 


2008.07.30 


05:20 


Deneb 


33.3 


10.3 


2008.07.30 


07:42 


Deneb 


32.3 


-9.9 


2008.07.30 


11:32 


Deneb 


30.1 


-37.6 


2009.07.27" 


07:51 


CALP 


31.8 


-19 


2009.07.27" 


08:11 


Deneb 


33.3 


-11.1 


2009.07.27" 


08:28 


CALP 


31.4 


-23 


2009.08.26 


06:51 


CALl" 


31.1 


-26.8 


2009.08.26 


07:06 


Deneb 


33.0 


-18.5 


2009.08.26 


07:25 


CALP 


30.5 


-31.0 



.10.01 


03:56 


Deneb 


33.2 


-14.3 


.10.01 


04:34 


CALl" 


33.2 


-14.3 


.10.01 


05:56 


Deneb 


32.2 


-27.3 


.10.01 


11:35 


Rigel 


25.3 


-6.5 



2009.10.25" 


02:21 


Deneb 


32.3 


-26.2 


2009.11.17 
2009.11.17 


01:35 
09:48 


Deneb 
Rigel 


32.8 
27.1 


-20.4 
-20.5 



" CALl: HD 184006, A5V, V=3.7, estimated diameter 

0.6±0.05mas from SearchCal@JMMC 
" Observations performed remotely from France 



vations were caiTied out at regular intervals between July 2008 
and November 2009. Rigel observations were performed in 
October-November 2009. The VEGA control system can be 
handled from France remotel y. A detailed description of this 
control system is presented in lClaussd (120081) . 

The red detector was centered around 656nm, and we 
made use of the high-resolution mode (R = 30 000). For each 
observation, VEGA recombined the SI and S2 telescopes 
forming the S1S2 baseline. This is the smallest baseline of 
the array, for which the extended Ha emission is not over- 
resolved. The S1S2 baseline is almost aligned north-south 
and the projection of the baseline onto the sky does not vary 
much during the night. This is considered an advantage in the 
context of our observations, which are designed to detect the 
time variability of the Ha line-forming region. Details of the 
observations can be found in Table |2] On several occasions, 
more than one acquisition of Deneb was performed during 
the night, providing interesting, albeit limited information 
about the spatial asymmetry of the source. Since our goal was 
to investigate the spectral and spatial properties of the Ha 
line relative to the continuum, emphasis was not placed on 
obtaining very precise calibration of the absolute visibilities. 
The angular diameters of Rigel and Deneb are well known 
(see Sect. 13. Il l, so we relied on the published values to scale 
the continuum absolute visibilities. Moreover, it is difficult 
to find a suitable bright calibrator when observing at high 
spectral resolution, which is not well-suited to an accurate 
determination of the absolute visibility which implies a large 
continuum window. Improving upon the current determination 
of the angular diameter in the visible would require a dedicated 
observing strategy using the medium-resolution mode, larger 



baselines coupled with a fringe tracking performed in the 
near-IR. This possibility is foreseen in the future. Therefore, 
calibrator observations were not performed systematically, but 
only to ensure that the instrument was behaving well. 

The data reduction method is fully described in 
( 120091) and we only summarize it shortly 



Mourard et al. 



here. The spectra are extracted using a classical scheme, of 
collapsing the 2D flux in one spectrum, wavelength calibration 
using a Th-Ar lamp, and normalization of the continuum by 
polynomial fitting. We note that the photon-counting cameras 
saturate when the rate of photons is too high locally. Because 
of the brightness of the sources, neutral density filters of 0.6 to 
1 magnitudes had to be applied, depending on the atmospheric 
conditions and the spread of the speckle images on the slit. 
Spectra with a signal-to-noise ratio (SNR) of 300-400 were 
routinely obtained, although clear signs of saturation were 
found in a few of them. After careful testing, we checked 
that the saturation had only a very limited effect on the 
interferometric measurements. When the quality of the spectra 
was such that telluric lines are observable, these are used to 
refine the wavelength calibration, although most of the time the 
VEGA spectral calibration was checked against the reference 
spectra obtained by amateur astronomers. A series of spectra 
for Deneb is shown in Fig.[T] 

The raw squared visibilities (V^) were estimated by com- 
puting the ratio of the high frequency energy to the low fre- 
quency energy of the averaged spectral density. The same treat- 
ment was applied to the calibrators, whose an gular diameter 
was computed using the software SearchCaU (^ Bonneau et al 



2006). The expected absolute visibilities of VEGA in the con- 
tinuum from the short baselines do not provide strong con- 
straints on the angular diameter and a calibrator was not sys- 
tematically observed. When such an observation is available, 
we carefully checked that the observed visibilities were com- 
patible with the expected ones. A contemporary measurement 
obtained with the medium resolution mode is presented in 
iMourard et alJ (l2009l) . 

The information in a line was extracted differentially by 
comparing the properties of the fringe between a reference 
channel centered on the continuum of the source, and a slid- 
ing scie nce narrow channe l, using the so-called cross-spectrum 
method (iBerio etal.ll 19991) . The absolute orientation of the dif- 
ferential phase was established by consid ering the change in the 
spectral slope of the dispersed fringes dKoechlin et alJll996h 
when changing the delay line position. Thus a positive value 
corresponds to a photocenter displacement along the S1S2 pro- 
jected baseline in the south direction. 

The width of the science channel was 0.02 nm in good at- 
mospheric conditions, and degraded to 0.08 nm in poor weather 
conditions. The rms of the spectrally dispersed visibility in the 
654-655 nm continuum ranges from 3-4% at visibility 1 in the 
highest quality nights (2008/07/28, 2008/07/30, 2009/10/01) to 
about 7-8% in medium nights (2009.07.27, 2009.08.26) and 
more than 10% in poor nights (2009/11/17). The rms of the 
diff'erential phases follows the same trend ranging from 1-2° 



http://www.jmmc.fr/searchcal/ 
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100 
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Fig. 1. Top: Montage of Ha observations at various epochs 
for Deneb and Rigel. We also include one observation of the 
calibrator (HD17643). The dashed lines indicate the spectra 
recorded by amateur astronomers. Bottom: Normalized multi- 
epoch Ha observations of Deneb. The color coding is the same 
as that used in the top panel. 



in good conditions, to 3-4° in medium conditions, and 5-6° in 
poor conditions. Figures. |2]-[3] illustrate the best observations 
of Deneb secured in 2008. A Gaussian fit was applied to the 
differential visibilities and phase to retrieve accurate informa- 



tion about the spectral FWHM and position of the interfero- 
metric signal. The spectral location of the differential visibil- 
ity and differential phase dips are stable at a level of 0.005 nm 
(~3 km s"'). Information from the blue camera was also used, 
as some important lines, e.g. Sin 6343-6371 A can be investi- 
gated (see Fig.[3]and Sect. 13.3b . 



2.2. Spectroscopy from amateur astronomers 

Several Ha spectra were obtained during the same period 
with the 0.28 m amateur telescope (Celestron 11) located in 
Castanet- Tolosan (France) equipped with the eShel spectro- 
graph and a QSI532 CCD camera (CCD KAF3200ME). These 
spectra were used in this study as an indication of the emission 
level and variability of the stars. The typical resolution of these 
spectra is -11 000. 

The reduction of these data was performed using the stan- 
dard echelle pipeline (Reshel software VI. 11). H20 telluric 
lines are removed by means of division by a synthetic H20 
spectrum using Vspec software - the telluric-line list was taken 
from GEISA database (LMD/CNRS). We coiTected for the 
diurnal and annual earth velocity are corrected for (spectral 
wavelengths are given in an heliocentric reference for a stan- 
dard atmosphere). Systematic differences were found in the 
VEGA/CHARA and amateur continuum correction producing 
the differences seen in Fig. [T] between the solid (VEGA) and 
dashed lines (eShel). 



3. Results 



3.1. The diameters of Deneb and Rigel 



For the most accurate estimates to date of the diameters 



of Rig el and Deneb, we refer to lAufdenberg et al.l (12008 , 
2006al) . in which CHARA/FLUOR observations in the K band 



with baselines reaching 300 m are described. These obser- 
vations infer a UD angular diameter of 2.76±0.01mas, and 
2.363+0.002 mas, for Rigel and Deneb, respectively. To first 
order, no significant change in these diameters inferred from 
the visible or the near-IR is expected, the wind being too tenu- 
ous to shift the continuum-formation region (see Sect.lU. This 
value of the Dene b UD angular diame t er agre es with the optical 
measurements of iMozurkewich et alj ( l2003h obtained with the 
Marklll interferometer using multiple observations with base- 
lines ranging from 3m to 28m with 2.34±0.05mas at 0.8 //m, 
2.26+0.06 mas at 0.55 yum, and 2.25+0.05 mas at 0.41 ^m. The 
few V^ measurements secured from the present observations 
with a very limited spatial frequency range agree with the 
Marklll measurements. We estimated the uniform disk diam- 
eter of Deneb using only the red camera to be 2.31+0.04mas 
at 0.65yum, and using both cameras found the diameter to be 
2.34±0.03mas. Only the best observations could be used, as 
the visual magnitude of the calibrator (V=3.7) was too faint to 
often obtain a reliable estimate of the absolute visibility, even 
using the broadest possible spectral band (i.e. 5nm). 
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Fig. 2. Top row: Normalized flux (upper curve) and dispersed visibility (lower curve; spectral band 0.02 nm) for the Ha observa- 
tions of Deneb obtained on 2008.07.28 (solid line) and on 2008.07.30 (dotted line). The resolution is 7? = 30000. Bottom: Same 
as top row, but now for the differential phase. A strong signal changing with baseline direction is observed, indicating that there 
was a significant asymmetry in the line-forming region at this time. 



Table 2. VEGA/CHARA V^ measurements performed on 
2009.07.27 and 2009.08.26. 



Date Wavelength Baseline 

[nm] [m] 



V' 



TP 

y- error 



2009.07.27 


654 


32.97 


0.421 


0.019 


2009.07.27 


657 


32.97 


0.444 


0.024 


2009.07.27 


630 


32.97 


0.376 


0.016 


2009.08.26 


654 


33.19 


0.418 


0.051 



3.2. The Ha line-formation region 

The Wa line is one of the most optically thick of all lines seen 
in the optical and near-IR spectra of hot stars, hence represents 
an excellent tracer of their winds. This line has mostly been ob- 
served by means of spectroscopy at various spectral resolution, 
and in some dedicated campaigns with a very intensive time 
coverage aimed at recoveri ng the t i me va riabi hty of the line- 
formin g region. We refer to iRzaevI (l2008h and Morrison et al 
(l2008h for the latest reports on Deneb and Rigel. This con- 
spicuous line-profile variability suggests that the wind itself, 
where the Ha line forms, is variable in its properties. For ex- 
ample, variations may take place in the ionization, the mor- 



phology, or the density structure of the wind. The significant 
changes in the line-profile shape with time indicates that the 
Ha line-formation region is asymmetric, and that this asymme- 
try changes with time. 

The Ha line observed in Deneb exhibited evidence of some 
activity during the 2008-2009 observations, but not of any dra- 
matic change in mass-loss rate. The Rigel spectra were also ob- 
tained during what appears to have been quiet periods anddiffer 



from t he "typical" spectrum shown in Fig. 1 of Morrison et al. 



(l2008h . As can be seen here in Fig.[Tl the Ha line is as deep as 
the carbon Cii lines at 657.8 and 658.3 nm in Rigel. This is not 
an unusual state, as Ha often appears in pure absorption and 
to be symmetric about the line center ( i n the rest-frame) about 
20% of the time (using lMorrison et al.1 (12008 !) statistics). 

Both the Deneb and Rigel dispersed visibilities exhibit a 
profound dip in the Ha line, suggesting that the line forms over 
an extended region above the continuum. These dips are sym- 
metric about the line center (in the rest-frame) in both objects. 
The FWHM of the visibility dip was estimated by perform- 
ing a Gaussian fit to data of the high quality nights in 2008. 
These measurements yield a FWHM of 0.215+0.007 nm, i.e. 
98+3 km s ' for Deneb (Fig.|2]i. By comparison, the visibiUty 



O. Chesneau et al.: The Ka line-formation region of Deneb and Rigel 



Radial velocity 
-200 200 400 600 






-200 


Radial velocity 
200 400 600 


LZ 




- - - -I- -,,,A~-V».W.~™~WWl,»,v^- ---- 


1.0 




\r—^: 


0.8 


- 


y 


0.6 


-A.^^ 


v?ywWWw4AA^': 


0.4 




Deneb : 


0.2 




- 


00 




2008.07.28.09.01 : 
PA=-10/-20 deg 



Radial velocity 
-200 200 400 600 




636.5 637.0 637.5 638.0 638.5 636.5 637.0 637.5 638.0 638.5 636.5 637.0 637.5 638.0 638.5 

30 L J 30 



20 



10 




20 



10 




-lo r I 1 -ioL_ 

636.5 637.0 637.5 638.0 638.5 636.5 
waveleneth in nm 



0.0 
63 
30 



2008.07.28.10.51 
PA=-37/-46 deg 



20 



10 




637.0 637.5 638.0 
wavelength in nm 



-lOL^ 
638.5 636.5 



637.0 637.5 638.0 
wavelength in nm 



638.5 



Fig. 3. Same as Fig. |2] but now showing Deneb observations obtained with the blue camera, whose spectral range covers the 
Sin 6371 line. For a clearer comparison, the ordinate ranges are kept identical. The Sin 6371 line is seen as a strong absorption 
(the continuum level of the spectrum is offset by 0. 1 for clarity), marginally resolved by VEGA/CHARA. However, no phase 
signal is observed. 



signal is narrower for Rigel, with a FWHM of 70+5 km s~'. 
One can estimate the visibility in the line by using a contin- 
uum derived from published values of the angular diameter. 
The corresponding UD estimates for the highest extension of 
the Ha line-forming regions are 4. 1+0.2 mas and 4.2+0.2 mas 
for Deneb and Rigel, representing 1 .75 and 1 .5/?*, respectively. 

The differential phases are a direct indication of the position 
of the Ha line-forming regions at different radial velocities, by 
comparison to the continuum considered as the reference of 
phase. For the 2008 data, the Gaussian fits provide FWHMs 
that increase from the baselines oriented near +10° to -40° from 
6 1+2 km s"' to 84+4 km s"'. Sub-structures are also detected, 
with the highest peak near the zero velocity, and two satellites 
at about 40km s '. The peak level reached by the differential 
phases follows a trend from a large photo-center shift at RA. 
close to 5-10° that decreases toward -30-50°. 

The differential phases recorded for Deneb between 2008 
and 2009 are shown in Fig. |5] Dramatic changes are observed 
with periods of large phase signal alternating with periods con- 
taining no detectable signal. 

The "calm" periods may potentially provide important in- 
formation about the rotation of the star. This characteristic sig- 



nal can be observed in the rotating circumstellar e nvironment 



of Be - star disks (see De l aa et al. in preparation; iBerio et al. 



19991: IVakih et all Il998t ISted 119961). and also directly o n 



the photosphere of a rotating star (iLe Bouquin et al.l l2009h . 



Approaching and receding regions of a rotating star experience 
different Doppler shifts and are thus spectrally separated. In 
the sky, the natural consequence is that the line-absorbing re- 
gions are seen by the interferometer at one or the other side 
of the continuum, generating the well-known S-shape signal in 
the differential phases. This signal might have been detected 
in Rigel data of the 2009-10-01 (see Fig. [Hi. The rms of the 
phase of these data is 7.7°, and the maximum and minimum 
of the signal are at 26° and -17°, respectively, at the blue and 
red sides at ~15 km s ' from line center. The signal is kept at a 
similar level when extracted with a double binning of 0.04 nm, 
and the phase rms is decreased to 5.6°. This is on the order of 
the estimated v sin / of the source of about 36 km s"' . We were 
unable to detect a similar signal for Deneb probably because 
of the limited spectral resolution, insufficient to resolve its low 
estimated vsini of 20km s"'. As the signal is anti-symmetric, 
it cancels out if the spectral resolution is too low. One must also 
keep in mind that the quality of the data fluctuated and that the 
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Fig. 4. Top: Ha observations of Rigel on 2009.10.01 (upper 
curve), and the corresponding differential visibility, scaled to 
the expected absolute level of the continuum (lower curve). 
Bottom: Differential phase curve exhibiting a clear S-shape sig- 
nal, indicative of rotating circumstellar material. A vertical dot- 
ted line is shown in both panels to indicate the location of the 
Ha rest- wavelength determined by a Gaussian fitting of the vis- 
ibility curve. 
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Fig. 5. Time evolution of the diff'erential phases of Deneb. 
The strong phase signal observed in 2008 is indicative of a 
very asymmetric environment, off-centerd from the continuum 
source. No clear S-shape signal is observed, in agreement with 
the low vsin/ of Deneb of ~20 kms ', unresolved by the in- 
strument. 



binning could not be kept identical at all epochs. It is not im- 
possible that weak signals are blurred (such a signal might be 
visible the 2009/11/17; Fig.|5]l. Given the large extension and 
the activity observed in the dififerential phases in Ha, this line is 
probably not ideally suited to inferring the rotation of the star 
To first order, and for marginally resolved sources, the dif- 
ferential phases can be considered to linearly depend on the 
photo-center position of the emitting source. This is no longer 
the case when the source is significantly resolved as in our 
case, but we consider this as a rough estimate. Assuming that 
the line emission represent a fraction f ratio of typically 50% 
of the hght in the Ha line core, the following relation is used 
p - -((f)/360){A/B)(l/fra,io)- We find that the astrometric shifts 
induced by the inhomogeneous Ha emission reach 0.5 mas 
(0.2/?*), but are lower than 0.3 mas most of the time. The event 
observed in 2008 is remarkable as the full line appeared to be 
off'-center relative to the continuum. A similar event was ob- 
served on 2009/10/01 to have an opposite direction. 

3.3. The Siu6347-6371A line-formation regions 

The Sin 6371 line was observed simultaneously with the Ha 
line in the blue camera (FigO]) in 2008. The line profile shows 
a pure absorption, but is slightly asymmetric with an extended 
blue wing due to extended absorption in the wind regions. The 
VEGA/CHARA observations of Deneb obtained in 2008 indi- 
cate that the line-formation region of the Sin 6371 Une is more 



extended than the continuum forming region, with differential- 
visibility dips of 10%, 19%, and 15% for baselines in the range 
[5°: 10°], [-10°:-20°], and [-35°:-50°]. The rms of the dispersed 
visibilities is 5%. The line FWHM estimated from Gaussian 
fitting is 55+3 km s"' and the FWHM of the visibility dip 
is narrower with a FWHM decreasing from 33+4 km s ' to 
24+4kms-i for R A. from [5°: 10°] to [-35°:-50°], respectively. 
One can roughly estimate the extension of the line-forming 
region using the minimum of the visibility and UD approx- 
imation to be 2.6 mas, 2.7 and 2.75 mas using the value of 
2.34 mas from NPOI as the diameter at 630 nm. One may spec- 
ulate whether this trend in the visibilities is permanent or a 
transient event closely related to the asymmetries inferred from 
the differential phases in the Ha line. Because of a spectro- 
graph realignment carried out in 2009, the Sin 637 1 line was 
not any longer observable with the blue camera in 2009, but 
the Sin 6347 was well-centered and could be analyzed. No such 
signal was observed at anytime. 

Despite the large depth of the silicon lines, no differential- 
phase signal was detected above an rms of <2°. This means that 
the imprint of the rotation of the star on the Sin line-formation 
region is not detected in the data, probably because of insuffi- 
cient spectral resolution. Strong Sin lines were also observed 
in Rigel's spectrum, but neither differential visibility nor any 
phase signal was detected. The SNR of the data is poorer than 
the Deneb observations, with an rms in the dispersed visibilities 
and phases of 10%. 
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Fig. 6. Left: In the bottom panel, we show a contour plot of the quantity P x I{P) as a function of Doppler velocity and impact pa- 
rameter P for Ho-, computed with CMFGEN using the stellar parameters of Deneb and a mass-loss rate value of 3 . 1 x 1 0^^ Mq yr" ' , 
corresponding to the best-fit model. This figure illustrates the distribution of intensity with impact parameter, and in particular 
serves to infer the regions that contribute significant flux to the line, and hence the spatial extent of the Ha line-formation 
region. In the top panel, we show the corresponding normalized synthetic flux in Ha. Right: Same as left, but now for a mass- 
loss rate value of 6.2xlO"^Moyr '. Notice the sizable change in profile morphology, echoing the change in the extent of the 
line-formation region. 



4. Comparison with radiative transfer models 



The BA supergiants, especially nearby ones, have been exten- 
sively studied with sophisticated radiative-transfer codes. It is 
not in the scope of this paper to determine more reliably the 
fundamental parameters of Deneb or Rigel. However, we wish 
to investigate several important issues related to the putative ef- 
fects of the wind on the interferometric observables. Deneb has 
been observed by several northern-hemisphere interferometers 
in both the visible and the near-IR. Rigel is currently moni- 
tored with VEGA/CHARA in the northern hemisphere, as well 
as in the southern hemisphere with the VLTI. Hence, in this 
study, we used the stellar parameters obtained by Przybilla et 
al. (2006) and Schiller & Przybifla (2008) for Rigel and Deneb, 
respectively. After developing a convergent model for a refer- 
ence mass-loss rate value at which Ha is predicted in absorp- 
tion, we then increased the mass-loss rate (keeping the other 
stellar parameters fixed) until the Ha line exhibited a well- 
developed P-Cygni profile. In this way, we explored the fol- 
lowings questions: 



- To what extent the angular diameter inferred from optical 
and near-IR continua is sensitive to changes in mass-loss 
rates? 

- Are the Ho-, Pa^, and Bry line-formation regions repro- 
duced by radiative-transfer simulations accounting for a 
wind? These lines form in different parts of the wind and 
can thus be used simultaneously to constrain its properties. 

In this paper, we focus on the Ha line because it can be ob- 
served with VEGA/CHARA. A similar study for Pa/? and Bry 
line-formation regions is postponed to another paper in prepa- 
ration. 

4.1. Numerical simulations 

Our radiative-transfer calculations were carried out with the 
Une-blanketed non-LTE model-atmosphere co de CMFGEN 
(lHillier&MiIleilll998t lOessart & HiUieill2005l) . which solves 
the radiation-transfer equation for expanding media in the co- 
moving frame, assuming spherical symmetry and steady-state, 
and under the constraints set by the radiative-equilibrium and 
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Fig. 7. Top: Theoretical visibility curves for Deneb, computed at selected wavelengths in the range ~6562-657 1 A, and thus con- 
taining the Ha line, and for mass-loss rate values of 3.1xlO"^Moyr"' (left) and6.2xlO"^Moyr"' (right). Bottom: Same as top, 
but now for Rigel and using mass-loss rate values of 1 (left) and 2x10"^ Mq yr ' (right). The dotted and dashed lines correspond 
to the best UD curves fitting the visible and the near-IR continua, respectively. The range of baselines of the VEGA/CHARA is 
indicated by two vertical dash-dotted lines. 



statistical-equiUbrium equations. It treats line and continuum 
processes, and regions of both small and high velocities (rel- 
ative to the thermal velocity of ions and electrons). Hence, it 
can solve the radiative-transfer problem for both O stars, in 
which the formation regions of the lines and continuum extend 
from the hydrostatic layers out to the supersonic regions of the 
wind, and Wolf-Rayet stars, in which lines and continuum both 
originate in regions of the wind that may have reached half its 
asymptotic velocity. 

We used the stellar parameters inferred by Schiller & 
Przy billa (2008) fo r Deneb, and by Przybilla et al. (2006) 
and iMarkova et al.l (|2008) for Rigel. The Mgii resonance 
lines suggest terminal wind speeds of ~240km s 'for Deneb 
(Schiller & Przybilla 2008) and -230 km s'for Rigel (Kaufer 
et al. 1996b). Their projected rotational velocities are low, i.e. 
20+2 km s 'and 36±9km s ' respectively. 

Even for sources as "close" as Deneb and Rigel, the dis- 
tance estimates remain quite inaccurate. Schiller & Przybilla 
(2008) derived a luminosity of 1.96x10^ Lq for Deneb, using a 
distance of 802±66 pc (assuming that Deneb is a member of the 
Cyg OB7 associatio n), while the one derived from Hipparcos 
dvan Leeuwenll2007h is significantly smaller (i.e., 432 +61 pc). 



which infers a luminosity estimate of 5.5+1x10"* Lq. This re- 
assessment directly scales down by almost a factor two the lin- 
ear scales of the size parameters. Furthermore, recalling that 
the mass-loss rate scales with the luminosity, this should imply 
a weaker steady-state wind mass loss for Deneb. The corre- 
sponding uncertainties are considerably lower for Rigel, which 
is much closer than Deneb. 

By adopting values these parameters from previous works, 
we find that the synthetic spectra computed by CMFGEN agree 
favorably with the observations. We thus adopt these stellar pa- 
rameters and vary the mass-loss rate value to assess the impact 
on the spectrum and in particular Ha and Sin 6347-637 1 .4A. 
In practice, we explore the effect of a wind mass-loss rate of 
1.55, 3.1, 6.2, and 12.4xlO"^M0yr-' for Deneb, and values 
of 1, 2, 4, and 8x10"^ Mgyr"' for Rigel. A typical eff'ect that 
appears in theoretical models is illustrated in the bottom pan- 
els of Fig. |6l where we show the distribution of the emergent 
intensity / (scaled by the impact parameter p), as a function 
of Doppler velocity and p. This typ e of illustration was in- 
troduced by iDessart & Hillien (l2005h to study line formation 
in hydrogen-rich core-collapse supernova ejecta. Here, it pro- 
vides a measure of the extent of the line-formation region of 
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Fig. 8. Left: In the upper panel, we compare the Hor line profile obtained in July 2008 (black; rectified) of Deneb and the 
corresponding model predictions for wind mass loss rate of 1 .5 (orange), 3.1 (red), 6.2 (green) and 12.4x 10"^ M© yr"' (blue). In 
the lower panel, we give the corresponding dispersed visibiUty for a nominal baseline of 33 m. Right: Same as before, but now 
for Sill 6371 A. 



He relative to the neighboring continuum and the sites where 
most of the emergent photons originate. Comparing the left and 
right panels suggests that a variation by a factor of two in the 
mass-loss rate value leads to sizable changes in the extent of 
the line-formation region, and correspondingly, large changes 
in the observed line profile. While spectroscopy is sensitive 
to the latter, interferometry is sensitive to the former Below, 
we describe the interferometric signals associated with these 
intensity maps computed with CMFGEN and produced using 
the stellar parameters suitable for Rigel or Deneb, and various 
mass-loss rates tuned to match observations. This is illustrated 
in Fig|7]in the case of Deneb and Rigel, using two different 
mass-loss rates. The visibility curves for various spectral chan- 
nels in the vicinity of Ha are plotted as a function of the spatial 
frequency. The differential visibilities observed by VEGA are 
generating by plotting the different value for a given projected 
basehne (between 27m and 33m). 



4.2. Visible and near-IR continuum forming regions 

In FiglTO] we show a zoom of the 0.6 yum and 2.2 //m squared- 
visibility curves of Deneb in the second lobe with the various 
mass-loss rate values used in this paper The visible contin- 
uum appears to be far more sensitive to such changes than 
the near-IR continuum. We note that these visibility curves 
are computed for a spectral resolution of 30 000, and can- 
not be directly compared with the broadba nd measurements 



of FL UOR/CHARA over the full /T'-band JAufdenberg et al. 
2008 h. 



Doubling the mass-loss rate from 3.1 to 6.2x10 ^Moyr ' 
does not significantly alter the optical thickness of the wind. 



The angular diameter of the star in the continuum near Ha de- 
termined by a UD fitting of the visibility curve does not there- 
fore change by more than 2%. However, we note that the sec- 
ond lobe is significantly affected. This conclusion is also true 
in the near-IR. All models also show that the near-IR UD an- 
gular diameter is systematically larger than the visible one by 
~2.5-3%, due to the increasing continuum opacity with wave- 
length (free-free and bound-free processes). The near-IR mod- 
els were scaled to match the accurate CHARA/FLUOR value 
of 2.363 mas, and this implies that the expected UD diameter in 
the Ha overlapping continuum is ~2.3 1+0.02 mas. This expec- 
tation agrees with the NPOI measurements of 2.34+0.05 mas 
at O.Sjum and 2.26±0.06mas at 0.55 ;um. 

For Rigel, the difference between the near-IR UD angu- 
lar diameter and its visible counterpart is increased slightly to 
reach -3.0-3.5%. Using the FLUOR/CHARA measurement of 
2.758 mas as reference, this would imply a diameter of about 
2.64 mas in the visible. In the near-IR, some instruments such 
as FLUOR/CHARA have an accuracy often better than 1% de- 
pending on the atmospheric conditions, and the ability to ob- 
serve with long baselines. The interpretation of the observa- 
tions obtained with this instrument might be affected by a mass- 
loss rate variation in the form of localized inhomogeneities, but 
probably on a smaller scale than an instrument with a similar 
accuracy operating in the visible.. 

The second lobe of the visibility curve is far more sen- 
sitive to any fluctuation of the mass-loss rate in the visible 
than the infrared. The reasons are twofold: a higher sensitiv- 
ity to limb-darkening effects in the visible, and a more ex- 
tended continuum-formation region, despite the very limited 
amount of flux involved (the wind remains in any case opti- 
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cally thin). This can be seen in FiglTO]as the second lobe of the 
visibiHty at 2.2/im is closer to the uniform disk model than the 
one at 0.6jjm. Balmer bound-free cross-sections increase from 
400 to 800nm. This causes the continuum photosphere to shift 
weakly in radius across this wavelength range and also alters 
the limb-darkening properties of the star. This eff ect might also 
be discernible in the Marklll data (IMozurkewi ch et al. 2003, 
Sect l3.1l ), Deneb appearing smaller at wavelengths close to the 
Balmer jump, although the baselines were too short to probe 
the second lobe of the visibility curve. An interferometer able 
to resolve a hot star up to the second lobe in the visible is 
very sensitive to the mass-loss rate, even in the case of a very 
weak wind. In the near-lR, the free-free emission strengthens. 
Free-free opacities increase at longer wavelengths, increasing 
the photospheric radius relative to that measured in the visible. 
Changing the mass-loss rate does not significantly affect the 
limb-darkening and the temperature scale near the star, thus 
the location of the extended emission. As a consequence, the 
free-free emission causes a larger diameter in the near-lR, in- 
dependently to first order of the mass-loss rate. 
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4.3. Dependence of the Ha and Sin lines on tiie 
mass-loss rate 

We then computed the spectrum and the dispersed visibili- 
ties in the Ha line at a spectral resolution R=30 000. We 
were unable to perform a satisfactory fit to the Ho- profile 
for either Deneb or Rigel, the absorption component being 
systematically too deep (Fig|8]l. This conclusion was reached 
in many studies, and these profiles could not be reproduced 
in terms of spherical ly-symmetric smooth wind models (see 
for ins tance Fig. 7 in|M arkova et al. ( 2008J), lAufdenber^ et al. 
(l2002h . and ISchiller & Przvbilla (.2008.) ). In this temperature 



regime, the models systematically predict profiles in absorp- 
tion partly filled-in by wind emission, hence only lower limits 
to the mass-loss rate can be derived by fitting the He line. 

However, this discrepancy is mitigated by the quality of the 
visibility fit shown in the upper panel of Fig|8] which is by far 
less sensitive to small absorption effects along the line-of-sight. 
Changing the mass-loss rate by a large factor of 2-4 has a dra- 
matic impact on the Ha line-formation region, changing the 
spectrum appearance and the dispersed visibility curve. When 
using a baseline in the range 26-33 m, there is a relationship be- 
tween the mass-loss rate and the FWHM of the visibility drop 
in the line that can be approximated by a second-order polyno- 
mial in the range considered (1.55 to 12.4xlO"^M0yr"'). One 
can see in Fig.fTTIthat the nominal model of Deneb with a mass- 
loss rate of 3. 1 x 10"^ M© yr"' fits the observed visibility curves 
well. The model FWHM of the visibility is 89 km s ' and the 
observed FWHM is 98+3 km s ' . Given the FWHM-mass-loss 
relationship, this would correspond to a mass-loss rate in the 
range 3.7+O.2xlO"^M0yr"'. For Rigel, a similar relationship 
provides a mass-loss rate in the range 1.5+0.4x10"' Moyr"'. 
We have studied the variations in 2008 and 2009 of the visi- 
bility profile without detecting any significant variation in the 
visibility profile, which, translated in a mass-loss rate variation. 



Fig. 9. Same as the upper panel of Fig. |8] but now for Rigel. 

suggests that the mass-loss rate has not changed by more than 
5%. 

By comparison, the Sin 6371 A line is less sensitive to any 
variation in the mass-loss rate (Fig[8]). The FWHM of the sim- 
ulated line weakly increases from 40 to 43 km s"', compared 
to the 2008 VEGA/CHARA measurement of 55±3km s"'. 
The FWHM of the coiTesponding model visibility dip are 
independent of the mass-loss rate, being unchanged at the 
value 19.5km s 'and lower than the measured mean value of 
29+5 km s"' for the 2008 observations. The depth of the visi- 
bility dip fits more accurately the curve with a mass-loss rate 
of 6.2 X 10"' rather than the nominal value of 3. 1 x 10"' Mq yr"' 
inferred from the Ha line. It is not possible to establish whether 
this issue is related to a particular event occurring during the 
2008 observations, or a permanent situation, or even a bias in 
the radiative-transfer model. 

Despite the lower quality of the interferometric data on 
Rigel, one can see in Fig|9]that the model provides a reason- 
able fit. We note that the Cii lines at 6578A and 6583A in the 
spectrum of Rigel are also slightly affected by the wind. Their 
visibility for the VEGA baselines are about 0.5-1% lower than 
the nearby continuum (using our model predictions), depend- 
ing on the mass-loss rate. Detecting such a weak signal would 
require that the accuracy on the differential visibility is about 
2-3 times better than the current instrument performance. 

4.4. A rough investigation of the morphology of the 
line-formation regions 

The variability in the differential phase signal observed is a 
clear sign of the stochastic activity of the wind of Deneb. Even 
though the data secured are partial, one can note that no dif- 
ferential phase signal exceeds 30°, which is impressively large 
considering that the baseline is limited. In 2008 and in October 
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Fig. 10. Synthetic continuum squared-visibility curves at 
0.6/zm (top) and 2.2//m (bottom) as a function of spatial fre- 
quency, and given for wind mass-loss rates of 1.5 (orange), 3.1 
(red), 6.2 (green), and 12.4x10"^ Mgyr^' (blue). The dotted 
black line corresponds to a uniform disk of 3.32mas, and the 
dashed black line to a uniform disk of 3.36mas. 

2009, the differential phase signal was able to be observed 
throughout the full line, at a level above 30°. This signal did 
not appear to be due to an increased level of stochastic 'noise' 
expected from the signature of small, localized clumps, but ex- 
hibited a well-structured signature. At other times, the differen- 
tial signal was only observed close to +vsim, and some patterns 
were reminiscent of the S-shape signature, originating in a ro- 
tating structure. 

The radiative transfer models demonstrate that at the in- 
ferred mass-loss rate of Deneb, significant emission originates 
in a circumstellar region of up to 2-4 stellar radii, and that the 
level of this extended emission depends strongly on the mass- 
loss rate. We tried to investigate the differential phases with 
ad hoc approaches using the models as a basis to generate the 
perturbation. Differential visibilities and phases provide con- 
tradictory and stringent constraints. On the one hand, a strong 
differential phase signal is observed, but on the other hand, the 
differential visibilities did not vary by more than 5% over the 
two years of sparse observations. This restricted considerably 
the size, flux, and location of the perturbation. The emergent 
flux of the best-fit model was perturbed using a 2D Gaussian as 
a weighting function. The parameters of the perturbation were 
its location, FWHM and flux ratio compared to the unperturbed 
model. The best parameter range was reached for a perturba- 
tion located at 2-3 stellar radii from the star, with a FWHM 
of 0.3-0.7 stellar radii, and a flux contrast of 10-25. One can 
see in Fig[l2] that the differential visibilities of the perturbed 
and unperturbed models are similar, whereas a significant dif- 
ferential phase signal is observed. The baseline is aligned in the 
direction of the perturbation observed in July 2008, i.e. north 
of the star. This model is also compatible with the differential 
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Fig. 11. Variation in the FWHM measured from the Yia vis- 
ibility curves computed with the intensity maps produced by 
CMGEN, shown here as a function of mass-loss rate and adopt- 
ing the stellar parameters suitable for Deneb. The adopted 
baseline is 33 m. For a high mass-loss rate value, the Ha 
line-formation region is fully resolved, and the FWHM of 
the visibility curve is close to the wind terminal velocity. 
The dotted lines indicate the mass-loss rate inferred from the 
VEGA/CHAR A measurements. 



phases observed with the two other baselines (Fig|2|l. Placing 
the perturbation farther away leads to a phase signal that is far 
weaker than to the observations, since the line-forming region 
of the unperturbed model naturally ends at 5-6 stellar radii. 
This is a limitation of this ad hoc approach. Making the per- 
turbation larger in size, even with a much lower flux contrast 
leads to a noticeable change in the differential visibilities. The 
same consequence is reached when putting the perturbation too 
close to the star. Finally, if the perturbation is too small (i.e. 
FWHM < 0. 1 mas), the flux contrast required to reach the ob- 
served phases is unrealistically high, and such that a significant 
'binary'-like perturbation in the visibilities should be observed. 
Finally, this perturbation study suggests that a large temporal 
variability in the visibility and phase is expected in the visible 
at high spatial frequencies corresponding to the second lobe of 
the visibility curve. 

Clarke & McLeanI (1 19761) and lMcLean & Clarkd (1 19791) did 



not detected any significant variation in the polarization of the 
HyS an d Ha lines in D eneb and Rigel, at a typical level of about 
0. l%. lHavesl (119861) performed an intensive broadband moni- 
toring in B-band polarization of Rigel, revealing a variability 
at a typical level of 0.2%. The variability pattern in the Q - U 
plane suggested that the ejection of material was not limited to 
a plane, and non-radial pulsation were thus proposed as the root 
cause of these localized ejections. This might be caused by the 
limited amount of intensive observations. This might also be 
a consequence, in the case of Deneb, of a small inclination of 
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Fig. 12. Comparison on the unperturbed model (dashed lines) with the perturbed model (blue lines, as explained in the text). The 
thick lines are the observation of July 2008. The parameters of the perturbation are the followings: flux contrast of 20, FWHM 
of 0.5 stellar radius, and position at 2.8 stellar radius south of the star. The insert shows the square root of the perturbed model at 
the wavelength corresponding to the core of the line, where the signal is maximum. 



the rotation axis on the sky. Any large-scale structures orbiting 
in the equatorial plane would generate a significant polariza- 
tion variability, without any preferential direction. In addition, 
no direction appears to be preferred by our observation, though 
we emphase that such a conclusion is very limited given the 
stringent limitations in the uv coverage. 



5. Deneb as a fast rotator 



Aufdenberg et alj (I2008ll2006ah presented evidence that Deneb 
is a fast rotator, based on 24 high accuracy FLUOR/CHARA 
visibilities in the /T'-band with projected baselines ranging 
from 106 m to 310 m, sampling the first and second lobes of 
the visibility curve. They detected slight departures from a 
purely spherical model at a level of about 2% in the near-lR, 
a discovery that had not been anticipated for an AB super- 
giant. They tentatively fitted the data with a rotating model at- 
mosphere, and demonstrated that despite the low v sin ; of the 
star, a model at half critical speed, seen nearly pole-on may 
account for the interferometric observations. A by-product of 
this fitting process is the determination of crucial parameters 
in this context, namely the inclination, estimated to be /=30°, 
and the orientation on the sky of the rotation axis found to be 
at about RA.=150° east from north. These findings have po- 
tentially important implications. In a general context, Deneb 
would be the first AB supergiant proven to be the descendant 
of a fast rotator. Support for this interpretation was provided by 



Schiller & Przvbillal (120081) . To more clearly interpret the pro- 
nounced mixing signature with nuclear-processed matter, these 
authors proposed that Deneb was probably a fast rotator ini- 
tially, and is currently evolving to the red-supergiant stage. 

In the context of our observations, the consequences on 
the Ha line-formation region must be evaluated. The Hot 
line may be affected by a moderate change in wind prop- 
erties, such as those that occur due to a latitudinal vari- 
ation in th e effective temperature o f the star, estimated to 
be -7 00 K ('Aufdenberg et al.! l2006ah . Moreover, as stated by 



Schill er & Przvbilla (2008), hydrogen lines are mainly sensi- 



tive to variations in logg. 

Neither the VEGA/CHARA data nor the theoretical study 
presented in this paper can provide definitive support for, or 
exclude, this interpretation. The Ha variabiUty observed by 
VEGA/CHARA is related to localized inhomogeneities in the 
wind of this star Is not impossible that that these inhomo- 
geneities may affect the continuum forming region in the K 
band, although is has been shown in Sect. I4.2l that large varia- 
tions in mass-loss rate are required to significantly affect the 
second lobe of the visibility curves. Yet, this conclusion is 
based on the ideal case of a spherical model, although given 
that the wind is very optically thin in the continuum, those in- 
homogeneities should not significantly alter the properties of 
the continuum-formation regions. 

The 2008 data arguably provide some support for the fast 
rotator model. We note that the S1S2 projected baselines were 
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roughly aligned (range of -45/+ 1 5 degrees) within the direction 
of the asymmetry found with CHARA/FLU OR at 150°, i.e. the 



directi on of the pole in the fast rotator model (lAufdenberg et al 



2006a[) . On the other hand, the Ha differential phases are ob- 
served to increase from a baseline roughly aligned to the pole 
direction at R\=150° . This may imply that the asymmetry is 
greater in a direction perpendicular to the pole. On the other 
hand, the Sin line dispersed visibilities seem deeper in the di- 
rection of the pole, suggesting a more extended line-formation 
region. Moreover, the high mass-loss rate inferred from the fit 
to the visibilities through the Sill line may also be an indication 
of a co-latitude dependence of the mass-loss rate. That no ev- 
idence of rotation was detected in the differential phases is an 
additionnal argument for a low v sin /, and therefore a nearly 
pole-on configuration for Deneb. These limited observations 
cannot provide definite conclusions, and the fast-rotator inter- 
pretation still needs to be investigated, both theoretically and 
observationally. This could be done, for instance, by repeating 
the FLUOR/CHARA observations to check whether the near- 
IR interferometric signal has changed or not, or by devoting a 
full VEGA/CHARA run with more extensive coverage, prefer- 
ably by using the 3T mode. It would also be of interest to in- 
vestigate the impact of the fast rotation of Deneb on the Ho- 
and Sin 6371 A line-forming regions, using a radiative transfer 
model of a rotating star with a wind. 

6. Conclusion 

We have presented the first high spatial and spectral observa- 
tions of two nearby supergiant stars Deneb and Rigel. The Ha 
line-formation regions were resolved and their angular size was 
found to be in agreement with an up-to-date radiative transfer 
model of these stars. The Ha line-forming region appears to be 
asymmetric and time variable, as expected from the numerous 
intensive spectral monitoring of AB supergiants reported in the 
literature. However, the time-monitoring of the dispersed vis- 
ibilities inferred from the Ha line of Deneb did not provide 
evidence of mass-loss rate changes above 5% of the mean rate, 
and the activity observed can be considered most of the time 
as a second-order perturbation of the wind characteristics. The 
increased coverage obtained in 2008 provides some evidence 
of a latitudinal dependence of both the Ha differential phases 
and the Sin 6371 line differential visibilities. This line was not 
significantly resolved subsequently. Only two observations of 
Rigel were secured. An S-shape signal was detected in the Ha 
differential phase of Rigel, suggesting that the rotation signal 
is detected. We note that the Ha line profile was almost pho- 
tospheric at the time of our observations. Observations with a 
large uv coverage may provide, in the case of Rigel, the direc- 
tion of the rotational axis on the sky and its inclination. 

Given the large angular size of Rigel and Deneb, the Ha 
line-forming region is fully resolved by the interferometer with 
baselines longer than 50m, and only the S1S2 baseline of the 
CHARA array is short enough to permit such an investigation. 
Therefore, it is not easy to significantly increase the observa- 
tional effort performed on these stars with such stringent ob- 
serving restrictions. An extension of this work is possible for 
sources with an apparent angular diameter between 0.5 and 



1 .5 mas, that are large enough, but also bright enough to make 
use of the highest spectral dispersion of the VEGA instrument 
(limiting magnitude of about 3). This concerns the AB super- 
giants closer than 1.5 kpc, and the supergiants in the Orion 
complex (d ~500pc) are in this context of particular interest. 
Simultaneous 3 telescope recombination would provide much 
better uv coverage than the one presented in this paper. Another 
interesting possibility is to add to this study several diagnos- 
tic lines such as the Can infrared triplet (849. 8nm, 854. 2nm, 
866. 2nm). These resonance lines are highly sensitive to non- 
LTE effects arising close to the photosphere and may shed some 
light on the regions were the material is launched. 
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